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The first ALMA view of IRAS 16293-2422 * 
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ABSTRACT 

Aims. We focus on the kinematical properties of a proto-binary to study the infall and rotation of gas towards its two protostellar components. 
Methods. We present ALMA Science Verification observations with high-spectral resolution of IRAS 16293-2422 at 220.2 GHz. The wealth of 
molecular lines in this source and the very high spectral resolution offered by ALMA allow us to study the gas kinematics with unprecedented 
detail. 

Results. We present the first detection of an inverse P-Cygni profile toward source B in the three brightest lines. The line profiles are fitted with a 
simple two-layer model to derive an infall rate of 4.5 x 10~ 5 M Q yr~' . This infall detection would rule-out the previously suggested possibility that 
source B is a T Tauri star. A position velocity diagram for source A shows evidence of rotation with an axis close to the line-of-sight. 

Key words. ISM: clouds - stars: formation - ISM: molecules - ISM: individual (IRAS 16293-2422) 



1. Introduction 

IRAS 16293-2422 (hereafter 116293) is a well-studied Class 
protostar with a bolometric luminosity of 32 L Q , embedded in a 
3 M Q envelope of size -3000 AU ( |Correia et al.|2004| >. It is lo- 
cated in the nearby p Ophiuchi star-forming region, at a distance 
of -120 pc pCnude & Hog|[T998l |Loinard et al.|[2008] i.I16293 
has been shown to consist of two main sources denoted as com- 
ponents A and B (hereafter 1 16293 A and I16293B), separated 
by 5" (600 AU) in the plane of the sky ( |Looney et aLl|2TjQ0) . 
The structure of I16293A might be more complex than that of 
I16293B: two centimeter sources (Al an d A2; |Wootten|[ l989 
and two submillimeter sources (Aa and Ab; Chandler et al. 2005 
were detected toward I16293A, while I16293B has not shown 
any sign of substructure at these wavelengths so far. 

Despite its low luminosity, 116293 has a rich chemistry, with 
hot-core-like (hot-corino) properties at scales of ~100 AU (e.g., 
Blak e et al.||1994HCeccarelli et al.|[T998l |Schoier et al |[2 002 ; 
Cazaux et al.||2003| |Caux et al.||2011t |Coutens et al.||2012V 
Single-dish observations and modeling (van Dishoeck et al. 
1993] [Ceccarelli et al.|[2000{ |Schoier et al.||2002) suggest that 
the emission from organic molecules arises from small-scale 
regions toward the continuum sources, where the temperature 
(~80-100 K) would allow grain-mantle evaporation. Moreover, 
recent interferometric observations have shown that some com- 
plex species are more abundant toward I16293A than toward 
I16293B ( [Bottinelli et al.|2004||Kuan et al. 



ward I16293B are much narrower (typically less than 2 kms~' 
wide). The determination of the centroid velocity toward each 
source is a matter of debate, where different studies present dif- 



ferent results. While Bisschop et al. (2008 1 argued that the high- 
excitation lines of the complex organi c s peak at Vj. sr =1.5-2.5 
km s _1 toward both sources, J0rgensen et al. ( 201 1 1 reported av- 
erage systemic velocities of 3.2 and 2.7 km s" L and average line 
widths of 2.6 and 1.9 kms~' for A and B, respectively. 



Several studies (e.g., Kuan et aL[ 2004; Bottinelli et al. 
Rao et al.||2009| |Bisschop et al.| 120081 ICrimier et al. 
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J0rge nsen et al.|201 1 p have focused on 116293 to characterize the 



respective nature of the two continuum sources. Indeed, 116293 
is known to have a quadrupolar outflow ([Walker et al.|| 1988 



Early interferometric studies (Bottinelli et al. 



2004). 



2004 Kuan 



et al.| 2004 > have shown that the spectra toward I16293A have 
broad lines (FWHMs up to 8 kms 1 ), whereas the lines to- 



* Continuum and spectral data are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http : //cdsweb . u- strasbg . fr/cgi-bin/qcat?J/A+A/ 



Miz uno et alJ[T990} |Castets et aLl|2TjOTj |Rao et al.||2009| l sug 
gesting the existence of a binary system. However, significant 
physical differences (chemical and kinematic) exist, suggesting 
different natures or different evolutionary stages of the two con- 
tinuum sources A and B. Whereas it is generally agreed that the 
1 16293 A component is protostellar in nature, it has been sug- 
gested that the I16293B component either represented a more 
evolved (T Tauri) star ( |Stark et al.|2004||Takakuwa et al.|2007| l 
or alternatively a very young object ( |Chandler et al.|2005p . 

Evidence for large-scale infall toward this protostellar sys- 
tem has been found using single-dish observations (Walker et al. 
1986[|Narayanan et al.|1998| l. |Chandler et al.| ( |2005) > suggested 
that there is evidence for infall toward I16293B from the ten- 
tative absorption feature seen in SO when imaged using only 
baselines longer than 55kA. 



2. Data 

116293 was observed on August 16-18, 2011 in ALMA 
band 6 as part of the ALMA Science Verification effort, 
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yielding a total on-source observing time of 5.4 hours. 
The array was used to point successively at two dif- 
ferent positions, (a,<5)j2 oo=(16:32:22.99,-24:28:36.100) and 
(16:32:22.71,-24:28:32.326), allowing the observation of two 
overlapping fields and obtaining an homogeneous response 
over the full extent of the binary system from the primary 
beam response of the ALMA 12-m antennas. The spectral 
setup consisted of observing the H2CCO ll(i,ii)-10(i,io) line at 
220.178 GHz, using one baseband centered on 220.182 GHz, 
with a channel width of 61 kHz and bandwidth of 234.375 MHz. 
Atmospheric variations at each antenna were monitored contin- 
uously using water vapor radiometers (WVRs), in addition to 
regular observations of a nearby phase reference source. Gain 
changes were tracked using regular hot/ambient load measure- 
ments 

The uncalibrated dataset was released, together with spec- 
trally averaged, low-resolution calibrated images, and is publicly 
available since April 20120 We calibrated and imaged the origi- 
nal interferometric data with the CASA0software version 3.3.0. 
Absolute flux calibration was performed using the CAS A Butler- 
JPL-Horizons 2010 model applied to Neptune observations, 
which results in an estimated flux uncertainty of ±10-15%. 
Bandpass calibration was performed using observations of the 
strong quasar J1924-292, while time-dependent gain calibration 
was derived by regularly observing (each 12 minutes) the nearby 
quasar J 1625-254. The calibrated measurement set was spec- 
trally binned to a channel spacing of 120 kHz (Hanning smooth- 
ing, 0. 16 km s _1 resolution), and then CLEANed using the Clark 
algorithm. The use of 16 antennas with a longest baseline of 
220-m results in maps with a synthesized beam size 2.2"xl.0"at 
220 GHz. The continuum image was produced by using line- 
free channels between 220. 103 and 220. 1 12 GHz, extracted from 
the visibilities table and imaged separately. After a first imag- 
ing of the continuum map, we performed self-calibration on the 
continuum data which lowered the rms noise in the continuum 
map by a factor of three. The resulting continuum map has a 
final rms noise level of 3.6 mJybeairT 1 with a 1.9"x0.9"beam 
(FigjTJ, while the theoretical rms in this map should have been 
~0.7 mJy beam 1 , if the noise had been completely decorrelated 
over the 175 channels used to build the continuum map. The con- 
tinuum was removed from the visibilities to produce the spectral 
cube, and we applied the self-calibration derived from the con- 
tinuum emission to the spectral data. The resulting rms in the 
line-free channels of the spectral cube is ~4.5 mJybeairT 1 per 
120 kHz-channel, while the theoretical noise is 4.3 mJy beam -1 . 

These ALMA data have a spectral coverage included in the 
SMA observations of |j0rgensen et al.| ( |2011| l, but the ALMA 
FWHM-beam is almost a factor 2 lower, while the spectral res- 
olution is increased by a factor 5 and the sensitivity by a factor 
15-20. 



3. Results 

3. 1 . Continuum emission and molecular complexity 

Sources I16293A and I16293B are clearly detected in the contin- 
uum map, shown in Fig. [T] The continuum peak fluxes of A and 
B agree with the fluxes reported previously by Bottinelli et al. 
( |2"054l at 230 GHz with a similar beam. The ALMA observa- 
tions also trace some extended continuum emission connecting 
sources A and B, which was not detected in the PdBI map of 



Botti nelli et al.| | |2"004) , and is most likely due to spatial filtering 
and poorer sensitivity of the PdBI observations. This extended 
continuum emission traces the inner region of the common pro- 
tostellar envelope as traced by single-dish observations, and in- 
cludes the position at which Remijan & Hollis (2006i detected 
an inverse P-Cygni profile in the low-density tracer lz CO (2-1). 
The continuum emission is used to estimate the mass as 
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1 http://almascience.eso.org/almadata/sciver/IRAS16293Band6/ 

2 http://casa.nrao.edu/ 



where we assumed optically thin emission and a dust opacity per 
dust mass (/C1.3 m J) of 0.86 cm 2 g" 1 (thick ice m antles coagulated 
at 10 5 cm -3 from |Ossenkopf & Henning|l994| and a gas-to-dust 
ratio of 100. The total flux, measured over the 10-cr region, for 
sources A and B is 1.30+0. 1 1 and 1.06+0.15 Jy, respectively, 
which implies a mass of 0.21 and 0.26 M Q for sources A and B, 
respectively. 

Thanks to the improved sensitivity and spectral resolution of- 
fered by ALMA, many (~50) molecular transitions are detected 
toward the two continuum sources, while previous SMA obser- 
vations at this frequency only detected six molecular lines to- 
ward each object (J0rgensen et al. 2011). The peak spectra to- 
ward source A and B are shown in Fig.[l] where we can appreci- 
ate the difference in the line profile toward these two objects. A 
thorough line identification and comparison of chemical proper- 
ties will be addressed in a future publication. 

3.2. Evidence for Infalling gas toward source B 

Toward I16293B, inverse P-Cygni profiles are unambiguously 
detected in three molecular lines (see Fig|2|and |C.l| l. These line 
profiles are clear evidence for infall in source B. Although an 
absorption feature in SO (76-66) toward source B was reported 
by |Chandler et al.| ( |2005| > as suggestive of infall, the ALMA ob- 
servations here analyzed are the first observations of the inverse 
P-Cygni profile toward this source. 

The inverse P-Cygni line profiles are modeled to extract the 
velocity information of the infalling gas and estimate the in- 
fall rates. Here we used the simple two-slab model described 
by Myers et al. ( 1996|l w ith the modification introduced by 
|Di Francesco et al.| ( |2001[ ) to take into account the continuum 
source. The model fits the infall velocity of the layers, V !n , its 
optical depth, to, velocity dispersion, o~ v , and excitation temper- 
ature of the layer on the rear, T r , while the excitation temperature 
of the foreground layer, Tf, is fixed at 3 K; details of the model- 
ing procedure are described in more detail in Appendix |A| 

The CH3OCHO-A molecular line with an inverse P-Cygni 
profile is shown in Fig. [2] with the best two-layer model shown 
in red (all three molecular lines with an inverse P-Cygni pro- 
file are shown in Fig. |C.l) . The blue excess in the line profiles 
is not fitted, and therefore the fit is mostly constrained by the 
absorption feature, see |Di Francesco et aTj ( |2001| i; |Kristensen] 
|et al.| ( |2012| l for similar procedures. Despite its simplicity, the 
two-layer model contains more parameters than can be fully con- 
strained because of degeneracies in the model, but allows a ro- 
bust determination of the velocity information. The parameters 
of the best two-layer model are listed in Table [T] 

The infall rate is estimated assuming spherical symmetry as 
M i„ fail = 47IT 2 H,„//mH Vm, where at radius r,-„ the infall velocity 
and density have values of V,„ and respectively; and [i is 
mean molecular weight of the gas (2.3). 

The infall radius, r,„, can be estimated assuming that the in- 
fall velocity is only free-fall, M = Vf n rj„/2G = 4nr? ni n /imn/3, 
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Fig. 1. Left: Continuum map of 116293 obtained at 220.1GHz with ALMA. The synthesized beam is shown in the bottom left 
corner. The rms noise is 3.6 mJybeairT 1 in the central 20" region of the map. Contours are drawn to -3, 3, 10, 22, 39, 61, 88, 120, 
and 157 time the rms noise, where negative contours are plotted using dotted lines. The peak continuum flux densities for sources A 
andB are 0.539±0.004 and 1.066±0.004 Jy beairT 1 , respectively. Right: Spectra of sources A and B are shown in the bottom and top 
panels, respectively. The rms noise in the line-free channels is 4.5 mJy beam 1 . The red lines (solid and dotted) show the molecular 



transitions previously identified by J0rgensen et al. (201 1 1, and the solid lines mark the transition where the inverse P-Cygni profile 
is found. 
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Fig. 2. Spectrum toward the continuum peak of I16293B for 
CH3OCHO-A. The red line shows the best two-layer model of 
infall fit. The best- fit model parameters are listed in Table [T] See 
Fig. C.l for fits of all three molecules. 



Table 1. Parameters for fitted "two-layer" model 



Line 


T, 


TO 


v in 


0-,, 




(K) 




(kms- 1 ) 


(kms- 1 ) 


CH3OCHO-E 


44±3 


0.48±0.04 


0.49±0.02 


0.30+0.02 


CH 3 OCHO-A 


46±3 


0.45±0.04 


0.49±0.02 


0.31±0.01 


H 2 CCO 


60±10 


0.33±0.05 


0.51±0.07 


0.39±0.03 



Notes. Vj„ is the infall velocity of the layers. t is the peak optical depth 
of each layer. <r v is the layer velocity dispersion. T, is the rear layer 
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Fig. 3. Intensity-weighted velocity for CH3OCHO-E. Contours 
are drawn to -3, 3, 10, 22, 39, 61, 88, 120, and 157 times the rms 
noise of the integrated intensity, 15mJybeam _1 kms _1 , where 
negative contours are plotted using dotted lines. The beam size 
is shown at the bottom left corner. See Fig. C.2 for velocity maps 
of all three molecules. 



The infall rates obtained using the best fits are 4.2, 
and 4.8xlO- 5 M yr 1 for CH3OCHO-E, CH 3 OCHO-A, 
H2CCO, respectively. 



4.5, 
and 



excitation temperature. Tt is the front layer excitation temperature. The Velocity m 3ps 



following parameters are kept fixed: Vu 
Tf = 3 K, and <b = 0.3. See Appendix[A|for a description. 



3.4kms-\ T c = 20 K, 



while also estimating the central mass, M, as that from a uni- 
form density sphere with radius r,„ and density n,„. Therefore, 
the accretion rate is estimated as 



3V 3 / V- \ 3 



,-t 



(1) 



The intensity-weighted velocity map for CH3OCHO-E is shown 
in Fig. [3] where the line-integrated intensity is also overlaid in 
contours (see Fig. |C.2| for maps of all three molecules). The mor- 
phology of the emission and its intensity-weighted velocity for 
CH3OCHO-A and CH3OCHO-E are quite similar, which is not 
surprising given the almost identical energy level and Einstein 
coefficients (see Table |B.2"] i. For H2CCO, which has a lower en- 
ergy level than CH3OCHO-A/E, the emission is more extended, 
but with a similar velocity structure. Source A presents a promi- 
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nent velocity gradient, with a velocity variation of ~5.5kms _1 
across the extent of the molecular emission, which is aligned 
along the axis connecting the submillimeter sources Aa and Ab 
previously detected by Chandler et al. (2005). In source B there 
is no evidence for such a strong velocity gradient, and in fact 
the velocity variation across the entire source B is ~0.4kms _1 , 
which is more than an order of magnitude lower than in source 
A. 

The position velocity (PV) diagrams for source A (along the 
direction shown in Fig. C.2 1 are consistent with rotation of a 
disk, which is close to being edge-on (to explain the strong ve- 
locity variation). The PV diagrams for the three molecular lines 
used are shown in Fig. C.3 This orientation would also agree 
with the orientation of the outflow seen in SiO by Rao et al. 
(2009). However, the present ALMA data do not rule-out the 
possibility of two velocity components at scales smaller than our 
synthesized beam, due for example to the emission of several 
sources in a multiple system (e.g. Aa, Ab; |Chandler et al.|2005) >. 

4. Discussion and conclusion 

The molecules where the P-Cygni profile is observed are com- 
plex organics, usually associated with the hot-corino and not 
with the ambient cloud (e.g. |Caselli et al.[T9 93) although some 
of them have also been detected toward outflows ( |Arce et aL] 
|2008| l. Therefore, we can rule out the absorption from the am- 
bient cloud or ambient envelope, and associate the observed P- 
Cygni profiles to infall from the inner envelope. This infall de- 
tection would rule-out the possibility that source B is a young 
T Tauri star, as previously suggested by |Stark et al. ( 2004[ l. We 
notice that it is only thanks to the improvement in spectral reso- 
lution and sensitivity of the presented ALMA observations com- 
pared to previous interferometric observations of this source that 
the inverse P-Cygni profile could be detected and modeled. 

The infall velocity, V,„, derived from the three different 
molecular lines is the same (0.50±0.01 kms -1 ) and supersonic. 
If the velocity dispersion of the absorbing layer is only caused 
by thermal motions (c s - 0.32 km s~' for gas at 30 K), the infall 
velocity is between 1.3 and 1.5x the sound speed, i.e. supersonic 
infall. The infall velocity and derived infall rates are consistent 
with the values derived by |Kristensen et aL ( 2012| l for a small 
sample of low-mass protostars. It is important to note that for 
source A we could not identify an inverse P-Cygni profile, if 
present, because of the steep velocity gradient. 

The PV diagram for source A is consistent with rotation seen 
almost edge-on and generated by a central object of O.53M , 
which is consistent with previous estimations using the source 
kinematics (e.g. Chand ler et al.|2003] l. However, this mass esti- 
mate is almost twice as high as the mass estimated using contin- 
uum emission (see also Rao et al. 2009). This discrepancy can 
be reconciled if the dust emission is partially optically thick at 
submillimeter wavelengths. 

The lack of a stronger velocity gradient in source B might 
be due to either rotation happening at smaller scales than those 
probed in these observations or because source B is closer to 
being face-on, which has previously been suggested based on 
the analysis of the continuum emission (Rodrig uez et al.| 2005; 
|Loinard et al."p007l >. 

This work illustrates the new capabilities of ALMA, which 
open a new window for studying the kinematics and chemistry 
in the early stages of star-formation. Using ALMA to carry out 
observations of the kinematics of source A at higher angular res- 
olution will for example allow one to distinguish between the 
two possible scenarios we propose here to explain the velocity 



gradient detected in the present observations, e.g edge-on disk 
or close binary system. 
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Appendix A: The two-layer model 



Table B.l. Position of sources 



Each inverse P-Cygni profile was fitted using the simple two- 
slab model described by Myers et al. ( 1996]> w ith the modifica- 



tion introduced by Di Francesco et a 



(|20bl|i to take into ac - 



count the continuum source, see also jKristense n et al.| ( |2012| l 
The model assumes that the central continuum source is opti- 
cally thick, emitting as a blackbody of temperature T c , and fill- 
ing a fraction of the beam, 3>. Two layers of gas, front and rear, 
are infalling toward the central source with an infall velocity and 
velocity dispersion V,„ and cr v , respectively. The rear layer is il- 
luminated by the background radiation, T/,. If for each layer the 
peak optical depth is to, then the expected line emission at ve- 
locity V can be expressed as 



AT B = [j f -J cr )[\-e- T >] 

+ {\-<S»(J r -J b )[\-e- {T >+ T ' ) \ , 

where 

/ cr = 07 c + (l-$)/ r , 
and 

~(V - (Vlsr + V in )f 



T f = t exp 
T r = t exp 



-(V - (Vlsr - V in )) 2 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



where Vlsr is the source velocity, and the radiation temperature 
is defined as 



Jx = 



[ex P (r /r,)-i] 



(A.5) 



where Tq = hvo/ks, and vo is the line rest frequency. 

Since the front absorbing layer is probably subthermally ex- 
cited, we assumed a conservative value for the excitation front 
layer of Tf — 3 K. The continuum temperature, T c , was chosen 
such that the continuum radiation temperature, J C (T C ), matches 
the peak continuum flux in the image. The continuum source 
beam filling fraction, O, which is not well constrained, was set 
to 0.3 to use the same value as |Di Francesco et al.| (|2001 ), who 
modeled interferometric observations with a similar angular res- 
olution as we did; this value is consistent with the line model- 



ing results of coarse angular resolution Herschel data Kristensen 



et al. (2012 1. Although the value of <I> used here is arbitrary, 
we have also performed line profile fits with higher values of 
<E> and the obtained infall velocity and velocity dispersion are 
unchanged, while the values for T r and To are increased and de- 
creased, respectively. Similarly, if a higher value of T c is used, 
T r and To are increased and decreased, respectively, while we 
obtain the same results for V,„ and cr v . Therefore our results are 
robust to the value of <t> and T c used. An initial set of models with 
the centroid velocity, Vlsr, as a free parameter were run, and all 
gave a Vlsr of 3.4+0.1 km s _1 . Given these results, we decided 
to fix the value of Vlsr at 3.4km s _1 and reduce the number of 
free variables in the fit. 

The x 1 was minimized using the mp fit fun procedure 
(Markwardt 2009 ), where a spectrum is generated for different 
parameters and then compared to the observed line profile. The 
following parameters were kept fixed during the minimization: 
Vlsr, T c , T f , and 0>. 



Source 


RA. (12000) 
(hh:mm:ss) 


Decl. (J2000) 
(dd:mm:ss) 


IRAS 16293 A 
IRAS 1 6293 B 


16:32:22.873 
16:32:22.623 


-24:28:36.50 
-24:28:32.49 



Table B.2. Spectral parameters for lines with inverted P-Cygni 
profile 



Molecule 


Transition 


Frequency" 


E h 


Ref. 






(GHz) 


(K) 




CH3OCHO-E 


17(4,13) - 16(4,12) 


220.166888 


103.15 


1 


CH3OCHO-A 


17(4,13) _ 16(4,12) 


220.190285 


103.14 


1 


H 2 CCO 


1 1(1,11) - 10(1,10) 


220.17757 


76.48 


2 



References. The spectroscopic references given are the most recent 
cited in the CD MS and JPL database s jPickett et al.|19 98 ; Miille r et af] 
[200T1[2005) . (l) |Maeda etaT1p008l l, (2) |Guarnieri & Huckauf|p003) 



Appendix B: Tables 



Table |B.l) lists the position of the continuum sources. Table |B. 2] 
lists the information for the spectral lines analyzed. 

Appendix C: Inverse P-Cygni profiles, Centroid 
velocity and position velocity maps 

In this appendix, we show the inverse P-Cygni profiles and their 
respective models in Fig. C.l The centroid velocity maps for 
all three molecules studied, CH3OCHO-A/E and H 2 CCO, are 
shown in Fig. C.2 The position velocity diagrams along the cut 
shown in Fig. C.2 are shown in Fig. C.3 
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Fig. C.l. Molecular line spectra from the ALMA data toward the continuum peak of I16293B. In each panel, a red line shows the 
best two-layer model of infall fit for each spectrum. The best-fit model parameters are listed in Table [T] 




16 h 32 m 23.5 s 23.0 s 22.5 s 
RA (J2000) 



Fig. C.2. Intensity-weighted velocity for the three lines studied, CH 3 OCHO-A (left panel) CH3OCHO-E (middle panel), and 
H2CCO (right panel). Contours are drawn to -3, 3, 10, 22, 39, 61, 88, 120, and 157 times the rms noise of the integrated inten- 
sity, 15mJybeam _1 kms _1 , where negative contours are plotted using dotted lines. The beam size is shown at the bottom left corner. 
The gray line is the cut used for the position velocity diagram shown in Fig. |C.3| 
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Fig. C.3. Position velocity maps of source A, along the direction shown in the rightmost panel of Fig. 



C.2 



for the three lines studied. 

From left to right: CH 3 OCHO-A, CH3OCHO-E, and H 2 CCO. Contours are drawn to -3, 3, 10, 22, 39, 61, 88, 120, and 157 times 
the rms noise, 4.5mJy beam 1 , where negative contours are plotted using dotted lines. The red dashed line and red solid circle show 
the 8 km s arcsec -1 velocity gradient and position of source A (Vz,s«=3.4km s ), respectively. Notice that in all panels an adjacent 
molecular transition line is marked with a vertical blue line, which can also be identified in Fig. [2] The spatial resolution is shown 
at the bottom right corner. 



